Near degeneracy and pseudo Jahn- Teller effects in mixed-valence ladders: 

the phase transition in NaV 2 5 
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We analyze the electronic structure of a mixed-valence ladder system. We find that structural 
anisotropy and complex electron correlations lead to on-rung charge localization and insulating 
character. Charge fluctuations within the rung of the ladder interact strongly to the lattice degrees 
of freedom, which gives rise to large pseudo Jahn-Teller effects. The phase transition in NaV^Os 
should be driven by this kind of mechanism. 
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Stoichiometric compounds in which the formal valence 
state of an element is fractional are referred to as mixed- 
valence (MV) systems. MV materials display remark- 
able properties, such as ferrimagnetism and a metal - 
semiconductor transition in magnetite, a transition from 
a ferromagnetic metal to an antiferromagnetic insulat- 
ing state in some half-doped manganites or the insula- 
tor-insulator transition towards a spin-gapped state in 
NaV^Os . The phase transitions in these compounds are 
often discussed in terms of charge ordering (CO) at the 
transition metal (TM) sites. In such models, electro- 
static interactions should yield at low temperature (LT) 
an ordered arrangement of the extra d electrons, those 
associated in the ionic picture with fractional oxidation 
states for the TM cations. For NaV20s it was initially 
believed that charge order is present at all temperatures, 
in the form of alternating chains of V 4+ and V 5+ 0. 
However, it was later found that above 34 K the vanadi- 
ums are all equivalent @>E| and some kind of CO should 
only occur in the LT phase, presumably with a zigzag 
arrangement of the d electrons 0, IU, [fj . In this Letter we 
investigate electron -lattice interactions in NaV^Os . We 
take explicitly into account the ligand degrees of freedom 
and find that most probably the phase transition at low 
temperature is induced by collective, pseudo Jahn-Teller 
(JT) effects involving TM-3d and oxygen- 2p hole states 
localized within V-O-V clusters. Our findings indicate 
different physics as compared to the commonly accepted 
d-d CO model. Similar JT type effects, we believe, may 
be relevant for other MV or doped TM materials. 

Corner-sharing [VO5] pyramids form two-leg ladders 
in NaV205 . Such ladders are assembled in quasi two- 
dimensional layers, with pyramids on adjacent ladders 
sharing edges of their bases. Interatomic interactions are 
strongly anisotropic, as shown by first principles band- 
structure calculations 0, S • This is related to obvi- 
ous structural features: V-O-V bond angles of approxi- 
mately 90° between adjacent ladders, shorter V-0 bonds 
across the rung than along the leg, and tilting of the 
[VO5] pyramids towards the rung bridging oxygen [3j. 
The insulating character is associated with this rung- 
leg anisotropy and relatively large inter-rung Coulomb 
interactions [3, H, Ifl S S S LLfll • In a first approximation 



a single d electron resides on each rung, in an essentially 
non-bonding d xy - d xy orbital, where the x axis is parallel 
to the rung and y is along the leg 0, ■ Alternatively, in 
a configuration-interaction like picture, a superposition 
of left and right occupied orbital states, dl xy d^ and 

^L,xy^k,xy > would correspond to each V-O-V rung. If 
the energy difference between the V 3d and O 2p levels is 
not too large, O 2p 5 configurations may contribute as well 
to the ground-state (GS) wave-function. The ligand-hole 
configurations with the largest weight would imply holes 
"localized" on the rung oxygen Or, d xy p y d xy , due to 
strongly anisotropic d xy -p x ^ y overlap integrals and lower 
Madelung potential at the Or, site as compared to the leg 
oxygens Ol 01 • Density functional calculations indicate 
an energy difference A = — e p sa 4 eV between the V d xy 
and 0-R.Py levels. However, neglecting inter-site Coulomb 
interactions, the on-site electron repulsion should reduce 
this number to A c s = A — U pp for p 6 d°^p 5 d 1 charge 
transfer. Due to the extended nature of the oxygen 2p 
functions there is no unambiguous procedure of estimat- 
ing the Coulomb repulsion U pp between two electrons in 
such an orbital. Nevertheless, semiquantitative estimates 
in copper oxides suggest a value as high as 4 eV 0] . 

In Table I we reproduce hopping matrix elements and 
orbital energies describing the d xy -p XtV and p x ,y~Px,y 
bonding on the V-0 ladder. These were obtained by ap- 
plying a powerful technique, the Vth-order muffin-tin or- 



TABLE I: In-plane NMTO orbital energies and hoppings for 
the high-temperature structure of NaV2 0s (eV). O l stands 
for an oxygen ion on the leg of an adjacent ladder. For the 



DMRG calculations we used e 
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Orbital energies 
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bital (NMTO) downfolding method recently developed in 
Stuttgart In order to generate an effective Hamilto- 
nian capable of modeling metal -ligand interactions, we 
downfold all basis functions other than V d xy and O 2p. 

We use the NMTO hoppings and orbital energies to 
investigate the nature of the many-electron GS wave- 
function of a V-0 Hubbard type ladder model. The 
Hamiltonian is written as 



<jj'>tr i 

+ Upp / I n^riji + V v< j n ia n ja/ . 

j <ij>aa' 

We neglect the apical ligands and restrict to a single or- 
bital at each site: V d xy , Ob,P v , 0^p x on the same ladder 
plus nearest neighbor (n.n.) 0' h p y orbitals on adjacent 
ladders. We also neglect the finite V d xy -OLP y over- 
lap, which should be zero for V and Ol atoms in the 
same zy plane, and 7r-type pp hoppings like those listed 
in the lower part of Table I . Still, we believe that our 
model is able to capture the essential properties of V-0 
ladders in NaV205 . GS expectation values are obtained 
by density-matrix renormalization-group (DMRG) meth- 
ods 0|- For these calculations up to 20 rungs on the 
V-0 ladder, open-end boundary conditions, and up to 
to = 2000 states to build the DMRG basis were used. 
The typical discarded weight was 10~ n — 10~ 10 . We em- 
ployed the code developed by Jeckelmann et al. . The 
results reported here were obtained with Udd — 3 eV, as 
deduced in 0,0, HI, an d the same U pp for Or and Ol- 

For reasonable values of the U pp and V p d parameters, 
the DMRG results show significant Or hole character for 
the GS wave- function, see Fig. 1(a). With U pp — AeN and 
Vpd = 0.25 eV, for example, the occupation of the rung p y 
orbital is 1.60, which indicates that the overall weights of 
the Or p y and p y configurations are nearly equal. Yet 
the ladder is insulating, with a finite charge transfer gap, 
see Fig. 1(b). The fact that the gap is strongly dependent 
on the rung - leg anisotropy was previously evidenced in 
|10|. Regarding the other parameters, the most impor- 
tant is U pp in our model. In simple words, for a given 
Vpd, the size of Udd "matters" when U pv is sufficiently 
high to give large occupation at the d sites. 

Physical intuition would point to a strong electron - 
lattice coupling for bond-stretching vibrations within the 
V-Or-V unit. The situation should be somewhat simi- 
lar to the in-plane O half-breathing mode in doped cop- 
per oxides, see for example ^(J and Refs. therein. To 
study such couplings in the MV ladder, we simulate the 
displacement of the rung ligands by introducinga mod- 
ulation of the rung transfer integrals, t r ~ r - 7 / 2 and 
inter-site Coulomb repulsion, V pd ~ 1/r. The effect of 
these displacements on the GS (electronic) energy is in- 
vestigated by DMRG calculations. However, before dis- 
cussing the DMRG results, we illustrate the nature of the 
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FIG. 1: (a) Occupation of the V d xy and Onp y orbitals. Up 
to the 2nd digit these numbers are constant along the ladder. 
For Ol Px and L Py the occupation numbers are not less than 
1.85. Results with L = 12 rungs and m=1600 DMRG states, 
(b) Charge gap A c for ladders with L = 12 rungs and 3L holes, 
where A c = l/2 [J3(3L+2) + E(3L-2) - 2E(3L)] . The values 
extrapolated to 1/L^0 through 1/L = 1/16 and l/L = l/20 
change by 15-25%, with larger differences for higher U pp . 



electron - lattice interactions on an oversimplified, single- 
rung model. We first assume that in NaV20s three elec- 
trons (holes) are "localized" on each d xy -p y -d xy "rung". 
For an isolated 3-site rung, we deduce then the expres- 
sion for the lowest energy eigenvalue corresponding to a 
superposition of p y a.ndp y configurations. With V p d — 0, 
A c g — A — Upp, and neglecting double occupation at the 
d sites, the Hamiltonian reads in matrix notation 
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where tL^ R (x) « i r /(l ±x) 7 / 2 , x — <Sx/r(VOR), Sx is the 
displacement of the Or ion along the rung direction, and 
t(VOr) w 1.83 A Qj. The lowest eigenstate is a doublet 
and has an energy E = 1 /2A cff - 1/2 ( A^ ff + tt\ + At R + 

4 (*i + t R- i l i fl) 1/2 ) 1/2 - We express E(x) as an one- 
dimensional (ID) series expansion and obtain: 



1 



E{x) = E(0) - - 



630 if 
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x 2 + 0(x A ) . 



We find thus a negative, "electronic" contribution to the 
harmonic force constant associated with this displace- 
ment. If we take the values from Table I, t r (0) =t r = 1.3 
and A = 3.5, choose U pp — 3eV, express the energy in 
units of meV and the distortion in units of percents of 
the VOr bond length, we find E(x)-E(0) w l/2fc cl x 2 , 
A; i = —11.7. This is a huge amount, about 35 eV/A 2 . For 
the isostructural CaV205 compound, a 3c? 1 Mott insula- 
tor [lg, the infrared active VOr bond-stretching phonon 
was associated to the B^ u mode at 515 cm -1 [19f. which 
gives a force constant k[ t — 15.6 eV/A 2 . For NaV20s , 
the appearance of continuum features in this frequency 
range, in both reflectivity B^, u and Raman B^g spectra, 
makes the interpretation difficult [lji|2ij]. Nevertheless, 
this continuum proves the existence of strong electron 
lattice interactions. 
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The large magnitude of this effect is confirmed by 
DMRG calculations on the V-0 Hubbard ladder, al- 
though more complex interactions are now involved. 
Contributions to the harmonic force constant induced by 
electron -lattice couplings are plotted for various U pp and 
V p d values in Fig. 2(a). The numbers were obtained by 
introducing an in-line pattern of (static) V- Or-V dis- 
tortions. Each Or, was displaced by the same amount 
and all other ions were kept fixed. To extract the k c \ pa- 
rameters we calculated E{x) for several small x values. 
We also considered the bending of the V-Or-V bond. 

Effects induced by the vibronic mixing of two or sev- 
eral electronic states are referred to as (pseudo) JT ef- 
fects [2l|. The presence of quasilocalized and relatively 
close in energy states plus the strong coupling between 
the electronic and nuclear motion in NaV20s can be in- 
cluded then within this class of phenomena. We checked 
for the existence of similar effects in CaV2 0s . CaV2 0s 
is also highly anisotropic However, the nature of 

the ground and low-lying excited states is different and 
also their energy separation is larger, not less than Udd ■ 
Such interactions are expected here to be weaker, which 
is indeed confirmed by DMRG calculations on a "half- 
filled" d 1 -p & -d 1 ladder. For U dd = 3eV 0, 2<U PP <5, 
and 0.25 < V p d < 0.75 the electronic energy per single 
rung varies as E{x) -E(0) « 1/2 k' cl 5x 2 , 2A< -k' el <7.0 
(eV/A 2 ). Although less drastic, this still should produce 
some phonon softening. We obtain thus a "bare" har- 
monic force constant k' = k' iT + \ k' cl | . We assume now 
that k 1 is also a reasonable estimate for the bare har- 
monic force constant in Nav^Os and display in Fig. 2(b) 
an overall force constant A = k' — \k c \ \ that should de- 
fine the 2nd order term in the ID Taylor expansion of 
the potential energy. We see that for certain values of 
U pp and Vpd A is negative, which shows that the system 
is instable with respect to such deformations. The finite 
amplitude of these displacements is related to the pres- 
ence of higher order, positive terms. In the simplest ap- 
proximation we arrive at the so-called </> 4 model, widely 
used to model local anharmonic effects in insulators [23| : 
V R = l/2Ax 2 + l/4Bx 4 . 

Near degeneracy and a large contribution of oxygen- 
hole configurations to the GS wave- function were actually 
predicted by multiconfiguration, quantum chemical cal- 
culations on small clusters [24|,|25j. Even more, these ab 
initio calculations pEj indicate a double- well potential for 
distortions of the rung oxygen along the x axis and give 
strong support for the model we propose here. Whereas 
the approximate treatment of so-called dynamic elec- 
tron correlation effects might somewhat overestimate the 
weight of the Or p y configuration and cannot guarantee a 
highly accurate value for the A parameter, we still can ex- 
tract from the results of [25j a good estimate for B, essen- 
tially determined by "hard" core -core repulsion. Using 
the units of Fig. 2(b) and Vg in - = V R {x ta 3%) = -42 K 
[H, we find B = 4Vg in -/x% = 0.18. The DMRG re- 
sults show a 4th order term in the purely electronic E(x) 
expansion of the order of 1% of this value. For the bond- 
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FIG. 2: (a) fe e ; (full lines, left scale) and k' el (dots, right scale) 
parameters describing the effect of rung distortions on the 
electronic energy per rung of the V-O ladder in NaV2 0s and 
CaV2C>5 (eV/A 2 ). (b) A parameter associated with the </> 4 
model, see text. Units of meV for energy and percents of the 
VOr bond length for the amplitude of the distortion were 
used, (c) V-V charge disproportion Snv, Sny = n(d,L,xy) — 
n(dR tXy ), for a zigzag arrangement of the rung oxygens with 
x = Sx/r (VOr) = ±3%. In each figure L = 12, m = 1600, and 
the same symbols were used to plot the curves. 



stretching vibration of the Or ion we can also evaluate 
inter-rung couplings. We model these oxygens by simple 
point charges and consider only Or- Or Coulomb inter- 
actions. To lowest order we find a n.n. interaction 



Vrr 



4 r 



3 = 1 



4) 2 



~r - ^ ) 



c\2 



k=l 



1 = 1 



where C b = g 2 r(VO R ) 2 /6 3 and C c = q 2 r(VO K ) 2 / c 3 are 
coupling constants within the be plane and a slightly more 
complicated expression corresponds to the inter-ladder 
coupling C a in the ab plane. With an "effective" charge 
for the Or ion q w 1.5 and the lattice constants from 
[3, these inter-rung couplings are C a — 0.19, Cb = 0.23, 
C c = 0.10. We note here that the DMRG GS energy 
for a zigzag configuration of the Or oxygens is always 
slightly lower than the energy for in-line displacements. 
This would give some extra contribution to the Cb pa- 
rameter. Parametrized Hamiltonians based on low-order 
Taylor expansions of the energy surface were intensively 
studied in connection with structural phase transitions 
in ferroelectrics, by both molecular-dynamics and Monte 
Carlo simulations |26l 1271 l2q . Whereas such a numeri- 
cal study is beyond the scope of this paper, we still can 
make a rough estimate of the transition temperature. In 
a mean field treatment T c ~ (2dC)A/B [2j|, where d is 
the system's dimension and C is an isotropic n.n. cou- 
pling. Although not accurate, this gives at least the order 
of magnitude of T c , as comparison to numerical results 
shows 2Q\. For U pp = 5 and V pd = 0.75, A= -0.51 and 
with C~ 0.18 we obtain T c m 47 K. We note, however, 
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that a number of other effects were neglected, such as 
charge polarization, other pd and pp hoppings, longer 
range electrostatic interactions, coupling to other modes, 
e. g. the V- Or-V bending, and coupling to the strain. 
The last effects should determine the actual LT crystal 
structure. 

We turn now to the analysis of the GS wave-function 
with displaced Or ions and plot in Fig. 2(c) the changes 
induced by rung distortions on the V d xy orbital occu- 
pation numbers. For the U pp and V p d values considered 
here, the on-rung V-V charge disproportion does not ex- 
ceed 0.31 for x — 3%. The values would probably de- 
crease somewhat when including on-rung polarization ef- 
fects. Recent interpretations of the experimental data 
pflj seem to agree with our results. 

In conclusion, we investigated the interplay between 
electron and lattice degrees of freedom for a MV ladder 
system. We found that structural anisotropy and com- 
plex Coulomb correlations lead to quasilocalized, three- 
hole states on each rung. Charge fluctuations within such 
an "unsaturated" metal 3d - O 2p- metal 3d structure are 
coupled to the lattice, which gives rise to strong JT type 
interactions. We suggest that the phase transition in 



NaV2C>5 is driven by cooperativity among these on-rung 
JT effects. We mapped results of DMRG calculations 
for a p-d Hubbard type ladder onto a 4 model. For 
reasonable values of the on-site and inter-site electron 
repulsion parameters we could reproduce in a mean field 
approximation the order of magnitude of the transition 
temperature. The spin-gap formation below 34 K is prob- 
ably related to the existence of different V- Or-V paths 
between consecutive rungs on the same ladder [g. l25ll30| . 
The dielectric anomalies at T c p| should be induced by 
zigzag ion ordering and only small V-V charge dispro- 
portion. We also mention that for certain values of the 
Coulomb repulsion parameters, see Fig. 1, the GS wave- 
function has significant oxygen- hole character. Strong 
contributions of oxygen-hole configurations were recently 
proposed for other MV oxides, such as the half-doped 
manganites, see [3lJ and Refs. therein. 
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